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Abstract 

The extension of the Boltzmann-Uehling-Uhlenbeck model of nucleus-nucleus collision is pre- 
sented. The isospin-dependent nucleon-nucleon cross sections are estimated using the proper 
volume extracted from the equation of state of the nuclear matter transformed into the form 
^ , of the Van der Waals equation of state. The results of such simulations demonstrate the de- 



pendence on symmetry energy which typically varies strongly from the results obtained using 
only the isospin-dependent mean-field. The evolution of the n/p multiplicity ratio with angle 
and kinetic energy, in combination with the elliptic flow of neutrons and protons, provides a 
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CN) | suitable set of observables for determination of the density dependence of the symmetry en 



ergy. The model thus provides an environment for testing of equations of state, used for various 
applications in nuclear physics and astrophysics. 
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Introduction 



One of the main goals of intermediate-energy heavy-ion collisions (HIC) is to study 
properties of nuclear matter, especially to determine the nuclear equation of state (EoS). 
HIC provide a unique possibility to compress nuclear matter to a hot and dense phase 
within a laboratory environment. The pressures that result from the high densities 
achieved during such collisions strongly influence the motion of ejected matter and are 
sensitive to the EoS. Within the hard work of the researchers over the three decades, 
the EoS of symmetric nuclear matter was studied in detail by the study of giant dipole 
resonances, collective flow as well as multifragmentation flfl. The EoS of isospin asym- 
metric nuclear matter is recently underway, particularly, for the density dependence of 
symmetry energy. Considerable progress has been made in determining the sub- and 
supra-saturation density behavior of the symmetry energy {5 -12]. The later part is still 
an unanswered question in spite of recent findings in term of neutron-proton elliptic flow 



ratio and difference 
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llj . However, the former one is understood to some extent [5- 8] . 



although, more efforts are needed for precise measurements. 

Transport model is very useful to treat heavy ion collision dynamics and obtain impor- 
tant information of nuclear matter EoS as well as the symmetry energy. In intermediate 
energy heavy ion collisions, the Boltzmann-Uehling-Uhlenbeck model is an extensively 



useful tool [13 



Q, 



which takes both Pauli blocking and mean field into consideration. 



The BUU equation reads 



d 3 p 2 d 3 p 3 dfl 



da 



NN 



dfl 



v 12 X [/,/ 4 (l - /)(1 - f 2 ) - ff 2 (l - / s )(l - U)} 



S 3 (P + P2 -Pa -Pa), (1) 

where f=f(r,p,t) is the phase-space distribution function. It is solved with the test 
particle method of Wong 15[ , with the collision term as introduced by Cugnon, Mizutani 



and Vandermeulen 16j . In Eq. ( CQ , 



dcrr 



dfl 



and i>i2 are in-medium nucleon-nucleon cross 



section and relative velocity for the colliding nucleons, respectively, and U is the single- 
particle mean field potential with the addition of the isospin-dependent symmetry energy 



term: 



U = ap + bp* + 2a s (— y Tz I, (2) 

Po 

where / = (p™ — p p ) / p, po is the normal nuclear matter density; p, p n , and p p are 
the nucleon, neutron and proton densities, respectively; t z assumes value 1 for neutron 
and -1 for proton, coefficients a, b and k represent properties of the symmetric nuclear 
matter while the last term, which describes the influence of the symmetry energy, can 
be obtained e.g. from simple Weizsacker formula, where a s represents the coefficient of 
the symmetry energy term and 7 is the exponent, describing the density dependence. 
Typical sets of mean field parameters cover substantial range, between the soft EoS with 
the compressibility K of 200 MeV (k = 7/6, ap = -356 MeV, bp* = 303 MeV), and the 
hard EoS with K of 380 MeV (k = 2, ap = -124 MeV, bp* = 70.5 MeV) [13]. 

It is the aim of the present work to estimate the effect of the symmetry energy 
parametrization within the equation of the state on the crucial component of the transport 
simulations, namely the in-medium nucleon-nucleon cross section. 



Isospin-dependent nucleon-nucleon cross sections 



When considering influence of the symmetry energy on emission rates of nucleons in 
nucleus-nucleus collisions, one needs to understand whether and how the medium rep- 
resented by the equation of state can influence relative probabilities of emission of pro- 



in-medium 
and 
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tons and neutrons. Theoretical investigations of the density-dependence o 
nucleon-nucleon cross section were carried out for symmetric nuclear matter 
significant influence of nuclear density on resulting in-medium cross sections was observed 
in their density, angular and energy dependencies. Using momentum-dependent interac- 
tion, ratios of in-medium to free nucleon-nucleon cross sections were evaluated via reduced 
nucleonic masses JiJ and used for transport simulations. Still, transport simulation are 
mostly performed using parametrizations of the free nucleon-nucleon cross sections, even- 
tually scaling them down empirically or using simple prescriptions for density-dependence 
of the scaling factor 20[ . In the present work, a prescription for estimation of the density- 
dependence of the in-medium nucleon-nucleon cross sections corresponding to the specific 
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form of phenomenological nuclear equation of state will be presented. Such possibility to 
establish a simple dependence of nucleon-nucleon cross sections on density, temperature 
and symmetry energy is potentially important for a wide range of problems in nuclear 
physics and astrophysics. 



Equation of state of nucleonic matter 

Based on the single-particle potential, shown in Eq. (J2]), one can construct correspond- 
ing equation of the state. Change of the pressure in thermodynamical equation of state, 
which is also a measure of non-ideality of a neutron or a proton gas, can be evaluated as 

dU 

^Pnon— ideal T \ T=const \y) 

dV 

where U is the thermodynamic potential, V is the volume and T is the temperature. 
When evaluating the thermodynamic potential It as a sum of single-particle contributions 
of neutrons and protons, given by equation ([2]), one arrives to expression 

P = Xn(^r4pT + ap 2 + bnp^ + 2 7 a sPo (-) 1+ ^) 

j3/2(Z n ) p 

+ XpiT^T^/tr + "P 2 + b*p 1+K - 2 7 a sPo (-) 1+ ^). (4) 
h/2(z p ) po 

where x n = p n / p and x p = p p / p are neutron and proton concentrations and Y(z) 
is the factor, a fraction of the Fermi integrals f n (z), assuring that Fermi statistics is 
taken into account. The parameters z n = p n /T, z p = p p /T are the values of fugacity 
of neutrons and protons, with p n , p p being the neutron and proton chemical potential, 
respectively. This expression appears to provide separate terms for pressure of neutrons 
and protons, which however can be combined to obtain the typical quadratic dependence 
on isospin asymmetry i". The resulting pressure is the weighted average between two 
terms, which can be, in similar manner to equation fl2]), summarily expressed as 

P = (^rl)pT + ap 2 + b K p 1+K + 2 7 a sPo (— ) 1+7 r 2 /. (5) 

f 3/2(2) po 

These terms can be interpreted as the equation of state of the system of particles with 
the corresponding single-particle potential, given by the equation (J2J). 
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Proper volume in the Van der Waals equation of state 



In order to find relation between equation of state and nucleon emission rates one can 
turn attention specifically to the Van der Waals equation of state. It can be written, using 
particle density p, 

( p + a'p*)(l-pb>) = (f^\)pT, (6) 

/a/2 (z) 

where the parameter a' is related to attractive interaction among particles and b' 
represents the proper volume of the constituent particles. In geometrical picture the 
proper volume of the particle can be directly related to its cross section for interaction 
with other particles. It is possible to formally transform the above equations of state 
for neutrons and protons (JS} (and practically any other equation of state) into the form 
analogous to the Van der Waals equation. Then, by comparison, one obtains values of 
coefficients 

a' = -a, (7) 

and 



_ bnp K + 2 7 a s (f o rr z I _ bap* + 2ja a (£)iT z I 



where the latter provides a measure of the proper volume of the constituents, nucleons 
in this case, as a measure of deviation from the behavior of the ideal gas. The proper 
volume of nucleon can be used to estimate its cross section within the nucleonic medium 

- = (^) 1/3 b ' 2/3 > (9) 
16 

which can be implemented into the collision term of the Boltzmann equation. 

Concerning the physical meaning of this procedure, for each point in the p-T plane 
the Van der Waals equation of state is found which behaves identically to the nuclear 
equation of state ([3]) in the vicinity of that point. Thus the dynamics of the system can 
be described using the two parameters of the Van der Waals EoS, of which one provides 



5 



a measure of the effective volume of the constituent at a given density and temperature. 
The variation of the constituent volume reflects the interplay of the long-range attractive 
interaction, leading to its apparent increase, with the short-range repulsive interaction, 
leading to its apparent decrease. 

The trend of the estimated values of nucleus-nucleus cross sections, obtained using the 
soft parametrization of symmetry energy with p 2 / 3 - dependence is shown in the Figure [1] 
for the temperature 20 MeV. The extracted values of the cross sections are surprisingly 
close to the expected value of geometric cross sections. With increasing isoscalar density 
the values of cross sections initially grow until they reach maximum values in the region 
around half of saturation density and then monotonously decrease. The increase at low 
densities possibly describes gradual deviation from the equation of state of ideal gas 
due to increasing attractive potential while at higher densities the decreasing trend may 
indeed represent the properties of short-range repulsive interaction. At high densities, 
with increasing isoscalar density the sensitivity to symmetry energy tends to decrease 
and around the density 2po it is practically lost, which however can be preserved using 
harder parametrizations of the symmetry energy. 

It is worth mentioning that a similar rise and fall of nucleon-nucleon cross sections 



was observed by the Aim, Ropke and Schmidt [18j and it was explained as a precursor 
effect of super-fluid phase-transition. Also in the present case this effect can be related 
to the phase transition in the nuclear matter, since it is caused by the same interplay 
of attractive and repulsive interaction which influences also the proper volume and thus 
extracted cross sections. 



Implementation into the Boltzmann equation 

While the Boltzmann equation (CQ) is formulated in terms of density, it does not explic- 
itly consider temperature. Therefore temperature T needs to be estimated independently. 
It is possible to estimate temperature using the Maxwellian momentum distribution of 
nucleons 

m = (2^T)37i e - T < (10) 
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FIG. 1: Extracted isospin dependencies of nucleon-nucleon cross sections for temperature 20 
MeV and various densities. 



where m is the nucleon mass. Using this formula, local temperature can be estimated from 
momentum distribution in the cm. frame by evaluating the momentum variance. At early 
stages of collision this can be done primarily for transverse momentum since it provides 
a measure of mutual thermalization of particles from the projectile and target, which 
proceeds by distant elastic collisions generating the transverse momentum. More violent 
collisions would lead to emission of colliding nucleons and thus would not contribute to 
thermalization of the source. This temperature estimate can be done without requiring 
stopping and formation of the source equilibrated in all three dimensions, closer analogue 
would be the friction of two dilute gas clouds passing through each other. 

Evolution of average temperature of the fireball with the time is shown in Fig. [2j The 

n n 

results were obtained using the Boltzmann-Uehling-Uhlenbeck equation (BUU) |13l . Il4j 
in reaction 48 Ca+ 48 Ca at 400 AMeV using the impact parameter 1 fm. Temperature was 
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FIG. 2: (Color online) Evolution of average temperature of the fireball (solid line) with the 
time. The total volume of the fireball is shown in arbitrary scale as dash-dotted line. The 
results were obtained using the BUU in reaction 48 Ca+ 48 Ca at 400 AMeV using the impact 
parameter 1 fm. Dotted straight line shows the fireball temperature estimate for a given beam 
energy obtained from the systematics of pre-equilibrium spectra 211 ]. 



evaluated for each time step in the cubic cells with the side of 1 fm. Average temperature 
was determined as a mean value of temperature over all cells where number of nucleons was 
sufficient (corresponding to density of p /10) and temperature thus could be evaluated. 
The total volume of this fireball is shown in Fig. |2] in arbitrary scale as dash-dotted 
line. Dotted line shows the estimate of fireball temperature for a given beam energy 



obtained from the systematics of pre-equilibrium spectra [21] . One can see that the 
average temperature over the fireball at its peak value exceeds the estimate from the 
systematics, while the value averaged over the lifetime of the hot fireball (between 10 - 
30 fm/c) appears to correspond to the value from the systematics. Thus it appears that 
the procedure introduced here leads to reasonable estimate of local temperature. 

Since the temperature, determined using the assumption of Maxwellian distribution 
represents the classical Boltzmann statistics, it can be corrected in order to reflect the 
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FIG. 3: (Color online) Comparison of the nucleon-nucleon cross sections in two variants of the 
BUU calculations. On the left panel are the isospin-dependent nucleon-nucleon cross sections, 
obtained as the proper volume of the Van der Waals form of the equation of state, as a function 
of density, while on the right panel are shown the corresponding nucleon-nucleon cross sections, 
obtained using standard energy dependent parametrization, used in BUU calculation. The 
results were obtained using the BUU in reaction 48 Ca+ 48 Ca at 400 AMeV using the impact 
parameter 1 fm. 

Fermi statistics, which fermions like nucleons obey. To achieve this, one needs to multiply 
the classical temperature T Boltz , corresponding to the Boltzmann statistics, by a factor 
( Tt \ )~ 1 an d thus the formula flHD will turn into 

b' 



'11) 



pTsoitz + + 27p a s (^) 1+7 T z /' 

which corresponds to classical case of Boltzmann statistics. Thus, remarkably, this 
classical expression can be used also for Fermionic (or even bosonic) particles, obeying 
their corresponding statistics. From practical point of view, in this way the non-trivial 
determination of the Fermionic temperature, depending of fugacity, can be avoided. 

Once the local temperature is determined, it is possible to implement the isospin- 
dependent nucleon-nucleon cross section, obtained using the formulas (TTTj) and fl9]), to the 
reaction simulation which can be used for determination of collision rate. The results for 
solution of such fully isospin-dependent version of Boltzmann-Uehling-Uhlenbeck equation 
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FIG. 4: (Color online) Relative difference of the isospin-dependent neutron-neutron and proton- 
proton cross sections as a function of isospin asymmetry of the volume cell. The results were 
obtained using the BUU in reaction 48 Ca+ 48 Ca at 400 AMeV using the impact parameter 1 
fm. 

will be presented in the next section. 
Reaction simulations 

The behavior of the in-medium nucleon-nucleon cross sections was investigated using 
the already mentioned reaction 48 Ca+ 48 Ca at 400 AMeV using the impact parameter 1 
fm. The soft equation of state was used, leading to incompressibility coefficient K =200 
MeV. For the isospin asymmetric part the "asystiff" parametrization was used with two 
symmetry energy terms, the kinetic term with the parameters a s i= 12.5 MeV and 71 = 
2/3, resulting from the Pauli principle, and the potential term with the parameters a S 2 = 
17.5 MeV and 72 = 2, respectively. 

Fig. [3] shows a comparison of the nucleon-nucleon cross sections in two variants of the 
Boitzmann-Uehling-Uhlenbeck simulations. On the left panel are shown, as a function of 
density, the isospin-dependent nucleon-nucleon cross sections, obtained using the proper 
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FIG. 5: (Color online) Evolution of the difference of n/p multiplicity ratio between VdW- 
BUU and fBUU calculations for three angular ranges. The results were obtained in reactions 
124 Sn+ 124 Sn and 112 Sn+ 112 Sn at 400 AMeV using the impact parameter 1 fm at the time 200 
fm/c. Squares show the result with stiff symmetry energy while triangles show results for soft 
symmetry energy. Soft nuclear equation of state is used. 

volume of the Van der Waals-like equation of state. On the right panel are shown the corre- 
sponding free nucleon-nucleon cross sections, obtained as energy dependent parametriza- 
tion of measured nucleon-nucleon cross sections 16] . It is apparent that while the isospin- 



dependent nucleon-nucleon cross sections essentially follow the 1/p 2 1/3 -dependence, the 
nucleon-nucleon cross section parametrization of Cugnon et al. leads to much larger 
spread, mostly due to its explicit energy dependence. Nevertheless, one observes that 
both parametrization cover essentially the same range of values of the nucleon-nucleon 
cross sections. Furthermore, from the comparison 22j of the parametrization of Cugnon 
et al. to in-medium cross sections at saturation density, calculated using the G-matrix 
theory by Cassing et al. (23 1 , it can be judged that the in-medium cross sections, ob- 
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FIG. 6: (Color online) Same as Figure [5] but for stiff nuclear equation of state. 

tained using the proper volume of the Van der Waals-like equation of state, are in better 
agreement with somewhat higher values of G-matrix in-medium cross sections of Cassing 
et al., which reflect properly the Fermionic nature of nucleons. 

In general, it is remarkable, that the in-medium nucleon-nucleon cross sections can 
be possibly directly related to the equation of state of the isospin asymmetric nuclear 
matter. This offers a more consistent description of the nuclear reactions and various 
astrophysical objects and processes in term of properties of nucleonic matter, expressed 
using the corresponding equation of state. However, one has to take into account that the 
equation of state of the isospin asymmetric nuclear matter describes isotropic medium 
and thus the extracted in-medium cross sections represent angle-averaged values. These 
values are used in this work, and compared to the results obtained with angle-averaged 



free cross sections of Cugnon et al. 



16[, and thus the effect of the equation of state 



of the isospin asymmetric nuclear matter on in-medium cross sections is demonstrated. 
However, one can consider possibility to implement angular dependence, either using 
the compatible microscopic calculations or from the observed experimental free nucleus- 
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nucleus cross sections. This possibility is beyond the scope of the present work and will 
be investigated in our future work. 
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FIG. 7: (Color online) Evolution of the difference of n/p multiplicity ratio between VdW- 
BUU and fBUU calculations for three angular ranges. The results were obtained in reactions 
124 Sn+ 124 Sn and 112 Sn+ 112 Sn at 400 AMeV using the impact parameter 6 fm at the time 200 
fm/c. Squares show the results with stiff symmetry energy while triangles show results for soft 
symmetry energy. Soft nuclear equation of state is used. 

The magnitude of the effect of isospin asymmetry on the nucleon-nucleon cross sections 
can be judged from Fig. HI which shows the relative difference of the isospin-dependent 
neutron-neutron and proton-proton cross sections as a function of isospin asymmetry of 
the volume cell. The results were again obtained using the BUU in reaction 48 Ca+ 48 Ca 
at 400 AMeV using the impact parameter 1 fm. One can see that the relative magnitude 
does not reach very high values even for most isospin-asymmetric cells, and the sensitivity 
of the Boltzmann-Uehling-Uhlenbeck simulation to the isospin-dependent nucleon-nucleon 
cross sections will result from the cumulative effect of large amount of nucleus-nucleus 
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collisions. 

The behavior observed in Fi gure s |3] and H] appears to be consistent with the behavior 
shown in Figure 4 of the work 19j. In both cases the absolute values of cross sections 
decrease monotonously with increasing density while the relative difference of the neutron- 
neutron and proton-proton cross sections increases with increasing asymmetry. Thus it 
appears that the procedure used to extract in-medium nucleon-nucleon cross sections by 
determining the parameters of corresponding Van der Waals EoS reflects the same physics, 
which is encoded even to a phenomenological equation of state such as Equation fl2]). 
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FIG. 8: (Color online) Same as Figure [7] but for stiff nuclear equation of state. 



Figure |5] shows evolution of the difference of n/p multiplicity ratios between Boltzmann- 
Uehling-Uhlenbeck simulation with both isospin-dependent mean-field and nucleon- 
nucleon cross sections which are correlated to each other by the Eq. [9] (based on the 
analogy with the van der Waals equation of state, thereafter we call this simulation as 
VdWBUU) and Boltzmann-Uehling-Uhlenbeck simulation with isospin-dependent mean- 
field and free nucleon-nucleon cross sections (thereafter we call it fBUU), for three angular 
ranges as a function of the kinetic energy in the center-of-mass system. The results were 
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FIG. 9: (Color online) Elliptic flow of neutrons and protons in reactions 124 Sn+ 124 Sn and 
112 Sn+ 112 Sn at 400 AMeV at the impact parameter 6 fm. Solid and open squares show results 
of VdWBUU calculation and fBUU calculation, respectively, with soft symmetry energy. Solid 
and open asterisks show analogous results with stiff symmetry energy. 

obtained using the BUU simulation in reactions 124 Sn+ 124 Sn and 112 Sn+ 112 Sn at 400 
AMeV using the impact parameter 1 fm at the stopping time 200 fm/c. Squares show the 
result with stiff symmetry energy parametrization (a a x— 12.5 MeV, 71 = 2/3, a S 2 = 17.5 
MeV and 72 = 2) while triangles show results for with soft symmetry energy parametriza- 
tion (a sl = 12.5 MeV, 71 = 2/3, a s2 = 17.5 MeV and 72 = 1/2). Soft nuclear equation 
of state (K = 200 MeV) is used in this case. Particles are considered as emitted when 
they are separated in the phase-space from any other particle and separation is large 
enough to assure that two particles are not part of a cluster ( a condition ApAr* > 2h 
is implemented ). One can see that implementation of isospin-dependent nucleon-nucleon 
cross sections leads to significant variation of n/p multiplicity ratio and this effect appears 
to evolve with both kinetic energy and polar angle. Variation of n/p multiplicity ratio 
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is more significant for the more neutron-rich system, which offers a strong argument for 
the use of neutron-rich exotic beams for studies of density dependence of the symmetry 
energy in the future. 

Figure IH] again shows evolution of the difference of n/p multiplicity ratio between Vd- 
WBUU calculation and fBUU calculation, in this case using stiff nuclear equation of state 
(K = 380 MeV). Also in the case of stiff nuclear equation of state one can see that 
implementation of isospin-dependent nucleon-nucleon cross sections leads to considerable 
variation of n/p multiplicity ratio, which again evolves with both kinetic energy and po- 
lar angle. Also here the variation of n/p multiplicity ratio is more significant for the 
more neutron-rich system. It thus appears that variation of n/p multiplicity ratio be- 
tween VdWBUU calculation and fBUU calculation provides a robust signal of the density 
dependence of nuclear symmetry energy. 

Figures [7] and [8] show results for difference of n/p multiplicity ratio between VdWBUU 
calculation and fBUU calculation analogous to figures [5] and El with the impact parameter 
set to be 6 fm. It can be seen that the effect of isospin-dependent nucleon-nucleon cross 
sections persists, with comparable magnitude, even in peripheral collisions. This offers 
possibility to study such signal of the density-dependence of the symmetry energy in a 
wide range of centralities and thus eventually to provide a strong signal of the density- 
dependence of the symmetry energy in a wider range of nuclear density. 

The fact, that the isospin-dependent nucleon-nucleon cross sections, obtained using the 
proper volume of the Van der Waals-like equation of state, need to be introduced in order 
to fully explore the isospin dependence in the Boltzmann-Uehling-Uhlenbeck simulations, 
directly affects the applicability of the symmetry energy parametrizations to the study of 
astrophysical objects such as neutron stars or supernovae. Due to increased sensitivity 
to isospin due to isospin-dependent nucleon-nucleon cross sections, the symmetry energy 
parametrizations may change significantly and that will affect the extrapolations toward 
the nuclear densities, typical for neutron stars and similar objects. On the other hand, 
increased sensitivity may offer more possibilities for the simulations of reactions of exotic 
nuclear beams, with the possible observation of stronger isospin-dependent signals. 

The recently performed simulations based on the Ultrarelativistic Quantum Molecular 
Dynamics model (UrQMD) suggest that one of the most promising probe of the strength 
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of the symmetry energy at supra-saturation densities is the difference of the neutron and 



proton (or hydrogen) elliptic flows 



ay 



25] . The simulations were performed using 



both stiff (7= 1.5) and soft (7 =0.5) symmetry energy parametrizations. An inversion 
of the relative strengths of the elliptic flow for neutrons and protons is observed when 
the symmetry energy parametrization is changed from the stiff behavior to the soft be- 
havior. Neutron and proton directed and elliptic flows were measured a decade ago in 
197 Au + 197 Au collisions at beam energies from 400 to 80 AMeV using the LAND neu- 
tron detector and the FOPI Phase 1 forward wall 
the UrQMD model |2J 
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271 ] . Comparison of predictions of 



provided a constraint on symmetry energy at supra-saturation 
densities from transverse momentum dependence of the neutron and hydrogen's elliptic 
flow parameter u 2 measured in the 197 Au + 197 Au system with FOP1+LAND, suggesting 
a value of 7 



0.9 ± 0.4 



10|, in agreement with findings at sub-saturation densities. 



However the statistics of these data set severe limits on the conclusions that can be drawn 
by comparison to transport model calculations. The promising results of re-analysis of 
FOPI+LAND experiment initiated proposal for a new experiment 29], which is one of 
the first dedicated explorations of the symmetry energy at high densities. The exper- 
iment uses the LAND calorimeter for neutron and charged particle detection, and the 
impact parameter is determined with a detection system with high effective granularity 
at forward angles consisting of several Csl rings of the CHIMERA multi-detector jsoj and 
the ALADIN Time-Of-Flight wall 31[ . In addition, flow of light fragments are measured 
with the Krakow telescope 32[ array positioned opposite from LAND. Uncertainties are 
expected to be reduced by a factor of 4-5, thus allowing to constrain theoretical calcula- 
tions. Expected results of this experiment will provide a welcome testing-ground for the 
model calculations introduced in the present work. 

It is of interest to estimate what effect will the introduction the isospin-dependent 
nucleon-nucleon cross sections, obtained as the proper volume of the Van der Waals-like 
equation of state, have on the resulting elliptic flow of neutrons and protons. Figure [9] 
shows calculated values of the elliptic flow of neutrons and protons (determined conven- 
tionally as second Fourier coefficient of the invariant triple differential distribution v 2 
relative to reaction plane) in reactions 124 Sn+ 124 Sn and 112 Sn+ 112 Sn at 400 AMeV at 
the impact parameter 6 fm. Solid and open squares show results of VdWBUU calculation 
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and fBUU calculation, respectively, with soft symmetry energy parametrization. Solid 
and open asterisks show analogous results with stiff symmetry energy parametrization. 
One can see that introduction of isospin-dependent nucleon-nucleon cross sections, vary- 
ing the ratio of neutron-neutron and proton-proton collision rates and thus modifying the 
in-plane/out-of-plane emission ratios, influences the resulting values of the elliptic flow, 
and since the effect appears to vary between neutrons and protons, it will strongly influ- 
ence the differential elliptic flow, thus making it a strong signature of the nuclear equation 
of state, as suggested in the work 10] . However, since the effect of symmetry energy on 
such differential elliptic flow tends to vary also with isoscalar part of the nuclear equation 
of state, it appears necessary to study differential elliptic flow in combination with other 
observables, such as the evolution of n/p multiplicity ratio at different polar angles and 
kinetic energies. Such a study, carried out on sufficiently neutron-rich system, can provide 
a good sensitivity to both isoscalar and isovector part of the nuclear equation of state. 



Conclusions 



The extension of the Boltzmann-Uehling-Uhlenbeck model of nucleus-nucleus collision 
is presented. The isospin-dependent nucleon-nucleon cross sections are estimated using 
the proper volume extracted from the equation of state of the nuclear matter transformed 
into the form of the Van der Waals equation of state. The results of such simulations 
demonstrate the dependence on symmetry energy which typically varies strongly from 
the results obtained using only the isospin-dependent mean-field. The evolution of the 
n/p multiplicity ratio with angle and kinetic energy, in combination with the elliptic flow 
of neutrons and protons, provides a suitable set of observables for determination of the 
density dependence of the symmetry energy. The model thus provides an environment 
for testing of equations of state, used for various applications in nuclear physics and 
astrophysics. 
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